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The mechanism of magnetization reversal in single-domain ferromagnetic particles is of interest in many 
applications, in most of which losses must be minimized. In cancer therapy by hyperthermia the opposite 
requirement prevails: the specific loss power should be maximized. Of the mechanisms of dissipation, here we 
0^ ' study the effect of Neel relaxation on magnetic nanoparticles unable to move or rotate and compare the losses in 

' linearly and circularly polarized field. We present exact analytical solutions of the Landau-Lifshitz equation as 

, derived from the Gilbert equation and use the calculated time-dependent magnetizations to find the energy loss 

per cycle. In frequencies lower than the Larmor frequency linear polarization is found to be the better source of 
^ ' heat power, at high frequencies (beyond the Larmor frequency) circular polarization is preferable. 
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^ ■ I. INTRODUCTION 

' The dynamics of the magnetization of single-domain ferromagnetic bodies under various conditions has been studied with 
. widely different motivations. Ferromagnetic resonance has provided perhaps the best known example. In this case a strong 
static magnetic field B generates and stabilizes the single-domain state. The torque G = /i x B exercised by this field, being 
perpendicular to B, causes a precession of the magnetic moment /i around the field, leaving the angle between /i and B intact. 
■ The resonance is detected by means of an alternating field perpendicular to B, whose energy is absorbed at a maximum rate at 
the Larmor frequency, ujl = 7|B|. 
r I Without a strong static field, in a sizable sample the single-domain state can be only achieved in a material of extremely large 
' O ■ permeability. This requirement dictated the choice of "soft iron" in the famous Einstein-de Haas experiment [1], which provided 
the first experimental determination of the gyromagnetic ratio 7. The limitation to soft magnetic materials can be circumvented 
though with very small samples. Clearly, samples smaller than the domain-wall thickness cannot contain more than one domain. 

The dynamics of the magnetic moment of nanometre-sized single-domain particles is of interest in connection with a number 
of applications. In magnetic recording, the processes taking place in magnetized particles under minor perturbations (reading) 
and during magnetization reversal (switching) are of paramount interest. In this case, the nanoparticles are immobile, their place 
and orientation are fixed. A fast process of magnetization reversal with minimal losses is the most important requirement for 
' applicability in magnetic recording, 
■r^lj- , Ferrofluids contain mobile magnetic nanoparticles, whose position and orientation can be controlled by a magnetic field. In 
T— I ' their application in hyperthermia, the most important requirement is the absorption of large losses under repeated magnetization 
I reversal. More precisely, the figure of merit to be maximized is the specific absorption rate (also called the specific loss rate, or 
heating rate), i.e., the amount of energy absorbed per second by a unit mass of magnetic nanoparticles^ The conditions under 
QQ which maximization has to be realized are quite demanding, as they are dictated by biophysical factors |2!]. 

. The mobility of magnetic particles in ferrofluids is an important feature in the context of losses under repeated magnetization 
^ ' reversal. As the particles are free to rotate, the potential energy E = — /i • B can be minimized in two ways: either the magnetic 
• ^ , moment rotates within the particle until it is aligned with the field, or the particle rotates as a whole. Two different processes are 
involved, Neel and Brown relaxation, respectively. In the latter the energy is transferred directly to the fluid, due to the viscous 
flow arising at the surface of the particle. In Neel relaxation the energy is transferred to the lattice, i.e., it is absorbed by phonons 
and magnons, before ordinary heat transport conveys it to the fluid. Both processes are characterized by appropriate relaxation 
times, ttv and tb, a parameterization that conceals the difference between the two mechanisms. The apparent similarity is 
further emphasized by the term 'magnetic viscosity', a misleading synonym for Neel relaxation. 

The common practice in hyperthermia is to expose magnetic nanoparticles present in a ferrofluid, preferably inside the malig- 
nant tumours, to a magnetic field alternating at a frequency of the order of 10^ Hz. The optimization of loss energy with respect 
to the amplitude and frequency of this linearly polarized field has been studied in detail However, the effect of circularly 
polarized field received relatively little attention. This neglect is not justified though, as there is no obvious reason why the 
specific absorption rate should be indifferent to the nature of polarization. 

In the present paper we report on the first stage of a systematic theoretical study of the dynamic response of a system of 
magnetic nanoparticles to linearly and circularly polarized magnetic field. The ultimate goal of this study is an analysis of 
the criteria for a maximum output of heating power applicable in hyperthermia. The relative importance of Neel and Brown 
relaxations is not clear. Hergst et al. [2] have made susceptibility measurements on a colloidal suspension of maghemite 
(7Fe203) particles at frequencies ranging from 10 Hz and 1 MHz. The very broad peak in the imaginary part of the susceptibility 
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found at 1 kHz disappeared when the particles were made immobile by letting the solvent freeze. In this case, viscous losses 
(Brownian relaxation) is clearly dominant. On the other hand, these authors also pointed out that the mobility of particles 
trapped in tumour tissue is not known. If they are clogged to the extent of not being able to rotate, the Neel mechanism is the 
only relaxation process available. Here we present analytic results for an immobilized single-domain particle. In Section 2 we 
briefly summarize the theoretical tools available for the study of the Neel mechanism and recapitulate the relation between the 
Landau-Lifshitz and Gilbert equations. Section 3 treats the case of linearly polarized magnetic fields, for which we present exact 
analytical solutions of the Landau-Lifshitz equation, which are used to determine the heating power as a function of frequency 
and amplitude of the alternating field. In Section 4 approximate analytic solutions are presented for circular polarization, which 
are exact under specific initial conditions. For arbitrary initial conditions, those solutions give an exact description of the 
dynamics in the steady state, which is stabilized after a short initial transient time-interval. The heating power is calculated for 
this steady state. 



II. NEEL RELAXATION 



We shall describe the magnetization processes in an isotropic single-domain particle by means of the Landau-Lifshitz equation 
101, which is generally used to consider the Neel relaxation mechanism ^ of magnetic nanoparticle systems. It is customary to 
write the Landau-Lifshitz equation in the form 

= /io7 f [m X H] + — [[m X H] X m] ) , (1) 
at \ ms J 

where m = n/V, V is the volume of the particle, /io = 47r x 10^^ N/A^ is the permeability of free space, 7 is an effective 
gyromagnetic ratio, a is a dimensionless damping constant, is the saturation magnetization and H = /Iq is the applied 
magnetic field. Formally, the Gilbert equation j^, 

= fiolQ ^[m X H] - 77[m x ^m]^ (2) 

is equivalent to the Landau-Lifshitz equation, as it can be verified |Q] with the identifications a = Ho^orjms and 7 = 7o/(l+a^) 
of the parameters of the latter in terms of those of the former Here 70 = 8.82 • 10^" Am^/Js is the gyromagnetic ratio of the 
electron and i] is the damping constant. As it was pointed out by Kikuchi HI, the factor 1/(1 + a^) appearing on the right-hand 
side of the Landau-Lifshitz equation as derived from the Gilbert equation removes a "physically implausible situation" 1^ in 
the predicted behaviour at large damping. 

The most conspicuous feature of the Landau-Lifshitz equation is that the vector standing on the right-hand side is perpendic- 
ular to the magnetization vector It follows then that dmn? jdt = m • dm/dt = 0, that is, the magnetization vector's magnitude 
does not change under the influence of the external field. Therefore, it is convenient to rewrite equation ([T]i in terms of the unit 
vector M = m/ms: 

^M = 7'[MxH]-a'[[MxH] xM]. (3) 

The new coefficients are given in terms of the ones in equations ([T]) and (|2]i as 7' = /io7 = /io7o/(l + ct^) and a' ~ fioja = 
fiQjQTjnis /(I + ct^)- With these definitions of the coefficients, equation (O is the Landau-Lifshitz equation as derived from the 
Gilbert equation 101, which will be referred to as the Landau-Lifshitz-Gilbert (LLG) equation in what follows. 



III. LINEARLY POLARIZED APPLIED FIELD 



In this section, we present an exact analytic solution for the LLG equation Q for an isotropic magnetic particle subjected to 
a linearly polarized, i.e. alternating applied magnetic field. The results enable the numerical evaluation of the heating rate per 
particle. The alternating field is assumed to be applied in the .T-direction, 

U ^ (cos(c^t),0,0), (4) 

7 

where ujl = 7'|H| is the Larmor frequency, oj is angular frequency of the alternating external magnetic field. 
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A. Free precession 



In the case of vanishing relaxation (a' = 0), the LLG equations for the Cartesian components of the magnetization vector are 

d 



dt 



0, 

cos{ujt) ujlMz, 
— cos{ujt) Lo^Ady. 

The analytic solution of these coupled equations for all possible initial conditions can be written as 



(5) 



My{t) = -y/l-M^o sin —sin{ujt) + 6a 



1 - M 



xO 



sm{Lot) + 5q 



(6) 



where M^i^ and are determined by the initial conditions: AI^o = Mx{Q) and 5^ = tan~^(Mj,(0)/M2(0)). The interpretation 
of equation (|6]l is that the unit vector representing the magnetization is precessing around the x axis with a time-dependent 
angular velocity. At time its projection on the yz plane is at an angle = ia,n~^(My{t)/Mz{t)) ~ [(o^i/w) sin((jjt) + (Jq] 
from the z axis. The angular velocity of the precession is changing periodically. If {lol/lo) >> 1, precession in a given sense 
prevails for many full cycles, before it slows down and changes sign. For the static limit, i.e., w — » 0, the solution (|6]l reduces to 

M^{t) = A-4o, 

My{t) = ^ 1 - A/Jo sin [ujLt + (5o] , 



Mz{t) = ^ 1 - AfJo cos [uLt + 5o\ . 



(7) 



This is the limit when the period of the alternating field is so long that the Larmor precession seems to proceed uninterrupted. In 
the opposite limit, [ujl/lo) « 1, there is only an oscillation of small amplitude about (Jq (i.e., &{t) k, (5o). 



B. Damped precession 



In the presence of Neel relaxation the LLG equations for the Cartesian components of the magnetization for linearly polarized 
applied field are 



d 
At 
d_ 
dt 
d_ 



A'4 = aNCOs{ujt) (1-MJ), 
My = cos(wi) (wlA^z - aNM^My) , 
= - cos{Lut) {cuLMy + aNM^M,) , 

where we have introduced ajv = Q;'|H| as a measure of damping. The analytic solution is 

(M,o - 1) + (M^o + 1) exp [2^ sin(c^t)] _ tanh sm{u;t)) + 



(1 - Mm) + [Mm + 1) exp sm{iut)] 1 + M^ tanh sin^)) ' 



My{t) = ^/l-M^(t) sin — sin(tjt) + (5o 



Mz{t) = ^yi-M^{t) cos —sin{uJt) + S^ 



OJ 



(8) 



(9) 



where, as before, M^ and Sq are determined by the initial conditions: M^ = ^/^.(O) and Sq = t&n^^ {My {0)/Mz{0)). The 
interpretation of equation (|9]l is that, like in absence of relaxation, the unit vector representing the magnetization is precessing 
around the x axis with a time-dependent angular velocity, but in this case the angle /? between the M vector and the x axis 
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changes as /? = cos~^{Mx{t)). At time t, the projection of M on the yz plane has a length of sin(/?) and is at an angle 
6{t) = ta.n^^ {My (t) /Mz{t)) = [{loi^/lo) sin(tjt) +<5o] from the z axis. The angular velocity of the precession and the projected 
length of the unit vector M are changing periodically. If (ajv/^) << land [ujl/uj) >> 1 the precessing unit vector approaches 
slowly the x axis, completing many cycles before it slows down and changes sign. Note that the first condition can be fulfilled 
both at very weak and very strong damping, since a^v oc (1 + a^)^^. For {un /i^) » 1 the precessing unit vector approaches 
quickly the direction of the driving field, it is very close to being aligned with the x axis, alternating between the positive and 
negative directions. For vanishingly weak or infinitely strong damping, that is in the ajv ^ limit, the solution reduces to 
equation (|6]l. For the static limit, i.e., lo Q the general solution (|9]l reduces to 

{Mm - 1) + {Mm + 1) cxp [2aNt] ^ tanh (ajyt) + M^ 

(1 - Mm) + {Mm + 1) cxp [2ajvi] ^ 1 + Mm tanh {uNt) ' 
v/1 sin [ujLt + So] , 

v/1 -Mj(t) cos [ojLt + So] . (10) 

In this case the M vector is precessing with a time-independent angular velocity ujl and converges towards the +x direction 
with a time constant ^ l/a^- In practice, for H w 10^ Am^^, we find that ~ jioloHa/ {1 + a^) w a/(l + cP') ■ 10^ s^^, 
so that even for weak damping, e.g., a = 0.001, the convergence is very fast. My and M^ vanish within a few microseconds. 
This implies that for frequencies much lower than 1 MHz we have M^ {t) = \. 



M,{t) = 

My{t) = 

MM) = 



C. Specific loss power 

The results given above enable us to calculate the energy loss for a single particle. The energy dissipated in a single cycle can 
be calculated as 

E = ,ol it[u.-y,o.^sl clt^H.— j. (11) 
Using the expression Q for H and the LLG equation dSjl and its solution Eq.®, the energy dissipated in a single cycle, 

E^l^iomsaNH dt [cos{ujt)]'^ [1 - M^{t)], (12) 

can be calculated numerically. The loss per cycle depends on the initial conditions, it has a maximum when AIm = 0. 

Analytic results are available in the limit of {um /^) << 1 where Mx{t) can be expanded in powers of QAr/tj. As in the 
initial state, in a ferrofluid, the magnetization vectors are randomly distributed, the relevant quantity will be the average of the 
loss power (fT2] i over this distribution. To lowest order in ajv/w this gives 



i-avcragc ~ -fJ-oT^msH . (13) 

O UJ 



At frequencies much lower than 1 MHz (aAr/tJ >> 1), the loss per cycle becomes independent of the initial conditions and 
frequency. In this limit it equals 

-^average ^ E ^ Adorns H. (14) 



IV. CIRCULARLY POLARIZED APPLIED FIELD 



Here, we present analytic solutions of the LLG equations for an immobile single-domain magnetic particle under circularly 
polarized, i.e. rotating applied magnetic field. Our exact analytic solutions are only valid for specific initial conditions. However, 
they also give the correct answer if one is interested in the behaviour in the steady state that sets in after a short relaxation process. 
The heating loss per cycle of a single particle in this state will be calculated to be compared in the next section to that obtained 
for the linearly polarized field. 

The appUed field is assumed to rotate in the xy plane with an angular frequency uj 

H = — (cos(wt),sin(cjt),0), (15) 

7 

We also introduce the angular velocity vector uj, which in this case is perpendicular the the xy plane. 
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A. Free precession 

The LLG equations for the Cartesian components of the magnetization are 

d , , 

— M:r = -u^lMz sm{ujt), 

= UJL Mz cos{ujt), 

= ujL{M^sm{ujt)-MyCos{Ljt)). (16) 

It is easy to find a special solution of these coupled equations for the case when d{Mz)/dt = 0, i.e., is a time-independent 
constant: 

M^{t) = cos{ujt), 

V ^ + ^1 

Myit) = ^^4==sinM), 

V ^ + 

Mz{t) = ^ (17) 

This solution is valid only for a special initial condition, viz., Mni;{0)/Mz{0) = lol/'-^, -^^a(O) = 0. For any other initial 
values, one has to find the general solution of equation ( fTSI l. This can be done by means of an appropriate rotation into a 
coordinate system in which the LLG equations have the simplest possible solution: a time-independent magnetization vector. 
The transformation is done in three stages. The first rotation, 

(+ cos(ti;t) + sa\{uot) \ 
-sin(wt) +cos(cji) (18) 
l) 

transforms the vector equation (fTSb into a coordinate system, which rotates around the z axis with the applied magnetic field. 
The transformed z axis points then in the direction of the angular velocity vector oj. The second rotation, 

-cos(e) -sin(e) 
O2 = I 1 I (19) 

+ sin(e) +cos(e) 

tilts the z axis into the direction of the angular velocity vector Vl ~ uj + uj^. As is aligned with H, it is perpendicular to lo, 
so that = = sjuj"^ + uj'j^ and cos(9) — lo/VL and sin(0) = lol/VI. Finally, the last rotation 

/ + cos{nt) - sin(m) \ 
03= +sin(m) +cos(m) (20) 

- V 01; 

is again a transformation into a rotating coordinate system. Similar rotations have been discussed in referencies ifToll and ifTlll . 
In the new rotating coordinate system the transformed magnetization vector (O3 • O2 • OiM) = (A^f, M^) is found to 

be time-independent. The inverse transformation (Oj^^ • 02^^ • 03^^) gives then the general solution of the LLG equation for 

circularly polarized applied field in the limit of vanishing Neel relaxation, which has the following form 



M^[t) = ^ cos(cjt) [M^ cos{nt) + M,, sin(f7t)] + sin(Ljt) [N^ sin(f7t) - cos{nt)] + cos{ujt)^l - A/| - M^, 
My{t) = ^ sin(tjt) [Mj cos(rjt) + M,, sin(m)] - cos(a't) [A/^ sm{Vlt) - M,, cos(m)] + ^ sin(wt).^l - - M^, 



Mz{t) = [M^ cos(m) + sin(f7i)] + '^^1- Aq - M^, (21) 

where A/^ and A/,, are the Cartesian components of the magnetization vector in the rotated coordinate system. They are deter- 
mined by the initial conditions as 



Af,= "^-°-"^V^^-J^■■^-°^^-^^^^")^ A/,.AV (22) 
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(M, and M^j are time-independent.) For non-vanishing relaxation, i.e. for ajv 7^ 0, = M^{t) and Mjj = M,i{t) become 
time-dependent. It is easily verified that = AI^ = implies the initial conditions for the validity of the special solution ( fTTj i 
and indeed equation (l2Tl i reduces to equation ( fTTI i in this case. In the limit of — ^ 0, equation (l2Tl i recovers the solution (|7]i 
which is obtained for static applied field with ajv ^ 0. 



B. Damped precession 



The LLG equations for the Cartesian components of the magnetization of a single magnetic nanoparticle in the presence of 
Neel relaxation, under circularly polarized applied field are 



-WL sm{ijjt) + UN [-M^My sm{ujt) + {My + M^) cos(wt)] , 

UL cos{ujt) + aM [-M^My cos{ujt) + + M|) sin(wt)] , 

^^L [Ma; sin(a;t) — My cos{u!t)] — [M^M^ cos{ujt) + MyM^ sin(cjt)] , 



(23) 



The general solution for arbitrary initial conditions can only be obtained numerically. It can be verified though that the LLG 
equations (|23]| are satisfied by the intuitively guessed solution 



Mx{t) = MxyC08{ujt - Lp) ^ M^^^""" cos.{ujt) 
My{t) ^ M.^y sa\{ujt -if)^ siii(wt) + M\ 



M^P^^siiiM), 



spec 



M.{t) = ^i-(A/:r)2 

where M^y is the projection of the magnetization vector onto the xy plane, M^q'^'^ 



M^y cos((p), M^l" 



(24) 

-M.jcy sm{if) and 



^4 



,2\2 



N 



2LoaN ^ ^ 

Equation (l24l i describes a precession around the z axis with a phase shift relative to the rotation of the magnetic field. However, 
there are conditions for the initial state: the phase shift is strictly determined by the tilt of the magnetization vector with respect to 

the z axis. The tilt is given by 61 = cos-^ ( ' ™^ *^ P*^^^^ ^^^^^ by ^ = tan"! (-M^p°7A/^p'''=). 
The relation between 9 and tp is thus implicit in equation (|25] |. 

We have verified with countless numerical calculations that any other solution of the LLG equation ( |23] | tends to equation 
(|24] | after a rather small transient time-interval. For the ensuing steady state the solutions of the LLG equation can be well 
approximated by equation (l24l i . Similar asymptotic solution of the general LLG equation for rotating applied field has been 
discussed recently by Sun and Wang ifioll . 

To understand the coupled dynamics of relaxation, precession and rotation it is useful to rewrite the LLG equation with 
Neel relaxation, equation ( l23T l, and the special solution ( |24] | in the rotating coordinate system defined in connection with the 
relaxation-free case. The transformed LLG equations read 







'-M,M,l 






'rM„M,^ 


5"< 




-M,M,^ 



-M„M5- cos(m) - M^M^-^ + (1 - M,;)- sin(m) 



n 

to 



-MqM^- cos(m) - M^Af^- sin(m) + (1 - M^)-^ 



(26) 



where (A/^, M^, M(^) represents the magnetization vector in the rotated coordinate system. As M remains a unit vector in its 
transformed form, the LLG equations for the components A/j, A/^ and Mq are not independent. The special solution (l24l) in 
terms of the rotating coordinates has the following form 



M^{t) 



M^jj cos{iU 



fa) 



Ml^""" cos(17t) 



-Af;rsm(^^t), 
M;P°^cos(m), 



(27) 
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where MIq" 



M^jj cos{ipa), -^,^0*"^ = ^M^Tj s'm{(pa) and similarly to equation {22\ . the time-independent parameters are 



, ,rSp0C 



(Miry 



, ,spoc 



, ,rSpOC 



(28) 



The advantage of the rotating coordinate system is that the time dependence that stands for the rotation is embedded inside the 
relaxation term, showing that the motion is dominated by the trend towards the stationary solution by a characteristic time con- 
stant ^1/ ajv- Therefore, it represents a good framework to look for the general solution of the LLG equation or, alternatively, 
one can try to find approximate solutions of the LLG equation. For example 



= M|P°^ sin(r!i) + M^P^^ cos{nt) + (i\/„o - AQ'") e 



s/2 . 



(29) 



well approximates the numerical solution of the LLG equation. The exponential decay of the last contributions to the in-plane 
components shows that the general solution converges quickly to the special solution. This is demonstrated in figure [T] for 
ajv = O.Ols^^, [a ~ 0.14) and uj = uji^ ^ 0.088 s^^, where the numerical and approximate solutions of the LLG equation are 
compared. The other consequence of the short transient process is, that irrespective of the initial conditions, the magnetization 
vector lies in the plane of the rotation, where it follows the driving field with a constant phase slip. 




200 300 

t[s] 



500 




200 300 

t[s] 



500 



FIG. 1 : Numerical (full line), approximate (dashed line), see Eq. l |29l l and the special (dashed-dotted line), see 
Eq. J27b solutions of the LLG equations for the rotated coordinates (t) and M( (t) of the magnetization are 
compared for a; = ljl ~ 0.088 s^^, a — 0.14. 



We note that both in the limits of weak and strong relaxation the phase slip ip goes to zero. This is seen in the values of the 
parameters in equation dZSl ) of the special solution ( l24b in the appropriate limits: 

M^rio^N ^ 0) ^ ^ , M;r(«A^-0) = 0, (30) 

and 

M:r(ajv ^ 108,s-i) « 1, M;r(«iv ^ 10«s-i) « 0, (31) 

as in the limit of strong relaxation on ^ 10^ s^^ values can be achieved. In general, the phase slip will depend on ujl{H) and 
u! according to the relation (p = tan^^ (— M^q°'^/A'/^q°'^). 



C. Specific loss power 



To calculate the energy loss per cycle for a single particle we use again equation ( fTTT i. The special solution (l24l i can be used 
to assess the merits of circular polarization, because in hyperthermia the field is applied for a long time (w lO'^ s) compared to 
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the duration of the initial transient process (« 10 ^ s). Using this special solution, the dot product in equation ( ITTT i turns out to 
be integrable and the energy loss per cycle is found to be 



(32) 



The weak-relaxation (or high-frequency) limit, << oj, ojl, can be studied by means of a Taylor expansion of (— A^^q"'^), in 
powers of + cj^). The result. 



E w 2irfj,QmsH 



(33) 



shows that there is no dissipation without relaxation. Making use of the relation apf — auuL this expression can be written in a 
form that highlights the scaling of the variables: 



E « 27rfj,QmsH — 



The loss per cycle for strong relaxation (or low-frequency), >> oj, ojl, shows similar scaling relations 



(34) 



E w 2TrfiomsH 



4 ? 



and 



E sa 2nfj,omsE[ — 



1 + [ay 



1 + 



(35) 



(36) 



We emphasize that the distinction between ujl and the 'bare' Larmor frequency uj^q — fiojoH cannot be ignored in the strong- 
relaxation limit. In fact, the definition = ojloC(/ (1 + a^) makes the limit aN ^ oo meaningless, because the largest value 
of ajy is uji^q/2, which is realized at a = 1. 



V. CONCLUSION AND SUMMARY 



We have reported numerical and analytical results on the response of an immobilised single-domain magnetic particle under 
linearly and circularly polarized magnetic fields. Comparison of the dynamics in the presence and absence of Neel relaxation 
enables an interpretation of the results in terms of an interplay between Larmor preccesion and relaxation towards energy minima. 




100k 1M 10M 100M 1G 100k 1M 10M 100M 1G 

Frequency (Hz) Frequency (Hz) 



FIG. 2: Energy loss per cycle (a) and specific loss power (b) for linearly and circularly polarised fields for 
magnetite as a function of driving field frequency for different a damping parameters. Arrows point towards 
the direction of decreasing a. Driving field amplitude H = 10'' A/m. 
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In figure 2a the resulting energy loss per cycle is depicted as a function of frequency and the dimensionless damping parameter 
a. Regarding the merits of linear or circular polarization of the magnetic field, we find opposite preferences for frequencies below 
and above the Larmor frequency. In the former case, which is relevant to hyperthermia, the energy absorbed by the magnet from 
linearly polarised field exceeds the absorption from rotating field by orders of magnitude. Apparently the relaxation towards 
a steady state with a small phase shift between the rotating H and m fields ensures a smooth process with very small losses. 
This cannot take place in a linear field, where the orientation of H changes abruptly twice in a cycle and relaxation towards the 
new energy minimum has to be repeated. In the case of high frequencies, which is not lacking technological relevance since the 
recent interest in materials with high microwave absorption II12I1 . the relaxation is slow compared to the rate of change of the 
magnetic field and the loss per cycle becomes frequency dependent. 

The tent-shaped curves in figure 2a reflect the relevance of the limiting cases worked out in sections 3.3 and 4.3. The low- 
frequency behaviour, which seems to prevail through many orders of magnitude in frequency, reflects equation ( fTOl i for the case 
of linear polarization and the first term in equation ( l34b for rotating field. Equation ( fTOl i describes a relaxation process, which 
takes place after each change of sign of H^, apparently at the same cost in energy as long as the period of the oscillation of 
the field is much larger than the relaxation time. The loss per cycle is then always the same and we see no dependence on the 
damping constant either, except that for weaker damping the deviation from the constant value begins at a lower frequency. 
Equation ( l34l i shows clearly that for frequencies much lower than the Larmor frequency the loss is always proportional to the 
frequency, the slope scaling with the dimensionless damping factor a (in the log-log plot this translates to straight lines with the 
same slope, shifted by distances scaling with a). In contrast, equations ( fT3] l and ( [33] l show that, in both cases, at high frequencies 
the loss is proportional to ajy /uj, which explains why all curves run parallel, with he same slope, and the dependence of their 
position on a is not monotonic. 

Figure 2b shows the same data in practical units, that is, it shows the proper specific absorption power in W/g for a single- 
domain particle of magnetite. At frequencies relevant to hyperthermia 10^ Hz) with linearly polarized field, respectable 
losses of kW/g order can be achieved, while the power absorbed from circularly polarized field is two to three orders of mag- 
nitude lower Our final conclusion is obvious: if Neel relaxation in isotropic sample is the dominant mechanism, the technical 
complications of generating a circularly polarised field in difficult geometry need not be considered. 
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